/* (1) Single outcome, multiple treatment groups */ data example; input treatment $ y; datalines; 1 102 1 104 1 101 1 105 2 105 2 106 2 106 2 107 3 112 3 108 3 110 3 109 ; run; proc glm data=example; class treatment; model y=treatment; lsmeans treatment/pdiff; lsmeans treatment/adjust=bon pdiff; lsmeans treatment/adjust=tukey pdiff; ods select diff; quit; proc glm data=example; class treatment; model y=treatment; lsmeans treatment/adjust=dunnett pdiff=control('3'); ods select lsmeans; quit; /* (2) Classical Multiplicity Adjustmet Procedures */ data two; input p @@; datalines; 0.01 0.03 0.20 0.52 ; proc multtest inpvalues(p)=two bon sid holm hommel hochberg; run;